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We compute spectra of particles produced during a dilaton-driven kinetic inflation phase within 
string cosmology models. The resulting spectra depend on the parameters of the model and on the 
type of particle and are quite varied, some increasing and some decreasing with frequency. We use 
an approximation scheme in which all spectra can be expressed in a nice symmetric form, perhaps 
hinting at a deeper symmetry of the underlying physics. Our results may serve as a starting point 
for detailed studies of relic abundances, dark matter candidates, and possible sources of large scale 
anisotropy. 
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I I. INTRODUCTION 

O 

. An inflationary scenario Jj],|| (the so called "pre-big-bang" scenario) postulates that the evolution of the Universe 
starts from a state of very small curvature and coupling and then undergoes a long phase of dilaton-driven kinetic 
inflation, in which the curvature and coupling grow, and then at some later time joins smoothly standard Friedman- 
Robertson- Walker (FRW) radiation dominated cosmological evolution, thus giving rise to a singularity free inflationary 
cosmology. Since the graceful exit transition from the dilaton-driven phase to the decelerated FRW evolution cannot 
■ occur while curvatures and coupling are small |}|4|] , an intermediate high curvature "string phase" j5j , in effect replacing 
O i| the big bang, has to separate the early inflationary phase from the late FRW phase. The required initial conditions 
of the Universe in this scenario were discussed ||-|^] , and the graceful exit transition is now better understood ||[l0) ■ 
We assume that the appropriate initial conditions were chosen such that a long dilaton-driven inflationary phase 
was indeed part of the evolution and that a later standard FRW cosmology followed, and wish to study possible 
consequences. 

We study production of a variety of particle types during the inflationary dilaton-driven phase. So far, most of the 
produced particle spectra were found to rise sharply with frequency. Therefore^in these spectra, most of the power is 
concentrated at higher frequencies, a property with interesting consequences , but also some disadvantages. In 

particular, at very large wavelengths there is almost no power, making these inhomogeneity perturbations an unlikely 
source for large scale anisotropy observed in the cosmic microwave background and required as seeds for structure 
formation. The standard explanation given for the generic spectral frequency dependence is that since the curvature 
increases sharply, particle production also increases and hence the resulting spectrum. Certain axionic perturbations 
were so far the only exception found to this tendency p4| . We find that although sharply rising spectra are indeed 
common, flatter spectra or even decreasing spectra are just as likely. We find that the slope of the spectrum depends 
on the spin of the particle and on the type of dilaton prefactor and whether it is massless or massive, revealing a rich 
range of spectral shapes, of which many more deserve further individual attention. 

The model of background evolution we adopt in this paper is a simplified model. The evolution of the Universe 
is divided into four distinct phases, the first phase is a long dilaton-driven inflationary phase. We assume that the 
background solution is the simplest solution of the string-dilaton-gravity equations of motion in 4 dimensions, the 
so called (+) branch vacuum. The second phase is a high-curvature string phase of otherwise unknown properties, 
followed by ordinary FRW radiation dominated (RD) evolution and then a standard FRW matter dominated (MD) 
evolution, both with a fixed dilaton. We assume that curvature stays high during the string phase, in a sense that 
we define better later, but we do not assume any specific form of background evolution. Obviously, by doing that 
we give up on the possibility of obtaining any information about the produced particles during the string phase. We 
postpone such analysis until some better ideas and more reliable methods for handling the string phase are at hand 
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(see, however, |^,|l5|]). We do not lose the ability to compute the spectra produced during the dilaton-driven phase, 
as we explain later. We parametrized our ignorance about the string phase background, as in 0, by the ratios of the 
scale factor and the string coupling, at moments in (conformal) time marking the beginning and end of the string 
phase zs — ai/as and gi/gs- It turns out that in the generalized setting we use, these two parameters are still 
sufficient to parametrize all spectra. 

Particles get produced during the dilaton-driven phase by the standard mechanism of amplification of vacuum 
fluctuations Jig , 17| . Deviations from homogeneity and isotropy are generated by quantum fluctuations around the 
homogeneous and isotropic background and then amplified by the accelerated expansion of the Universe. In practice, 
we compute particle production in the standard way, by solving linearized perturbation equations with vacuum 
fluctuations boundary conditions. The general solution of the perturbation equation is a linear combination of two 
modes, one which is approximately constant for wavelengths large compared with the curvature of the Universe (outside 
the horizon), and one which is generically time dependent outside the horizon. We understand the appearance of a 
constant mode as the freezing of the perturbation amplitude, since local physics is no longer active on such scales. 
The existence of the time dependent mode can be most easily understood in terms of a constant mode of the conjugate 
momentum of the perturbation |Is|| . The amplitude of the conjugate momentum also freezes outside the horizon, since 
local physics is no longer active on such scales. This forces the time derivative of the perturbation to be non- vanishing, 
leading to a "kinematical" time dependence of the perturbation amplitude. The analysis of perturbation equations 
and their the solutions in the Hamiltonian formalism was first suggested in Jl8| , as a way of understanding certain 
symmetries of the spectra of gravitons and photons. 

In addition to the dilaton and metric, string theory contains many other fields that have trivial expectation values 
in our model of background evolution and do not affect the classical solutions, but they do fluctuate. We are interested 
in "low-mass" fields, which are either massless, or have masses much below the string scale, such as moduli, gauge 
bosons, and their superpartners. In realistic models the mass of different fields may depend in a complicated way on 
the dilaton or other moduli. We will take the mass as an additional parameter and assume that it vanishes in the 
dilaton-driven phase and takes a constant value, much smaller than the string mass M s = M p g\ (M p is the Planck 
mass), from the start of the RD era and on. This assumption could be relaxed without substantial changes in results. 
For simplicity, we assume that the produced particles arc stable particles, decoupled from the radiation thermal bath. 
Adding the effects of interactions with the thermal bath and decay is rather straightforward and has to be done on a 
case by case basis. 

For a large class of particles, it is enough to specify a particle type with two parameters, m and I, specifying the 
interaction of the field with the background metric and dilaton. We explain the appearance of these parameters in 
more detail later, but for now we give the following list for some particular cases. For the dilaton or graviton m = 1 
and I — — 1, for the model independent axion m = I = 1, for perturbative gauge fields such as photons, m = and 
/ = — 1, while for some non-perturbative gauge fields m = and I = 0, Ramond-Ramond axions and non-perturbative 
scalars have m — 1 and 1 = 0. In general, fields with more tensor indices tend to have smaller (more negative) m's. 
We do not know at the moment a general rule for allowed values of m and I, but it is obvious that many combinations 
are possible. 

Possible benefits of our results are expected to come from several different aspects. First, from the different 
dependence on the fundamental parameters gs, z$ of decreasing as well as increasing spectra. This leads to different 
constraints on the allowed range of parameters, allowing to narrow their range and give better prediction for the 
spectrum of relic gravitational radiation background. Flatter spectra allow the possibility of producing large scale 
inhomogencities with enough power to explain observed cosmic microwave background anisotropics. Another exciting 
possibility is that some of the weakly interacting particles will constitute the dark matter in the Universe. Our results 
provide a starting point for actual calculations regarding these issues and we hope to study them in the future. 

Until now, the spectra that were computed were the following. Production of gravitons, for which the full dependence 
on parameters, effects of late time evolution and observability were taken into account p|,pT|. Production of dilatons 
was partially studied p3| , not taking into account the effects of an intermediate string phase. Photon production was 
studied fl2]| , taking into account dependence on parameters and late time evolution. Certain axionic perturbations 
were studied without taking into account the effects of an intermediate string phase and the effects of mass |l4j . In 
|l9f some additional fields were studied. 

Because this work is the first attempt to make a comprehensive survey of particle production in string cosmology 
models, in the hope of discovering some interesting phenomena, as we indeed do, we adopt a modest, minimalist 
approach. Our main goal in this paper is to show that spectra of produced particles come in different shapes, and 
we feel that this point is better demonstrated with spectra than can be reliably computed. We use approximations 
that we judge not to hide our main results and are accurate enough for our purposes. There are, however, many 
obvious improvements, which could be implemented if necessary. We would like to issue a warning about numerical 
coefficients, which should not be taken too seriously, in view of the approximations we used. In the same spirit we 
make no attempt to obtain bounds from constraints on the total energy in fluctuations, or to require that the total 
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number of particles leads to the observed radiation (as in (l9[| ), or to impose any other additional requirements on 
the spectra. We believe that this should be done for each class of particles separately, with additional input in the 
form of assumptions about stability, interactions and dependence on late stages of the evolution. We are completing 
a study of relic abundance axions f20|j , and plan to study some additional interesting cases. 

The plan of the paper is as follows. In section 2 we describe in more detail the model of background evolution, 
the perturbation equations, and their boundary conditions. In section 3 we describe our method of solving the 
perturbation equations, and solve them. The paper is quite technical and we have made a special effort to summarize 
our results in a self-contained form in section 4. The last section contains our conclusions and some preliminary 
consequences. 



II. BACKGROUND EVOLUTION AND PARTICLE PRODUCTION 



In this section we describe the model of background evolution, the perturbation equations obeyed by the different 
fields, and our method of solving these equations. 



A. The model of background evolution 



The model of background evolution we adopt in this paper is a simplified model. The evolution of the Universe 
is divided into four distinct phases, the first phase is a long dilaton-driven inflationary phase, the second phase is a 
high-curvature string phase of otherwise unknown properties, followed by ordinary FRW RD evolution and then a 
standard FRW MD evolution. We assume throughout an isotropic and homogeneous four dimensional flat Universe, 
described by a FRW metric with the line element ds 2 — a 2 (if) (drj 2 — 5ijdx % dx^ , where rj is conformal time and a(rj) 
is the scale factor. The dilaton <j> is time dependent, <fi = 4>(if). All other fields that we will be interested in are 
assumed to have a trivial VEVs. Therefore, to specify the phases in a concrete way, we need to specify two functions, 
a(rf) and <p(rj) and the conformal time boundaries between the phases, which we proceed to do. 

• The dilaton-driven inflationary phase lasts while i] < rj s , and while it lasts the scale factor and dilaton are given 
by the solution of the lowest order string-dilaton-gravity equations of motion, the so-called (+) branch vacuum, 

/„\(l-V3)/2 /n\-^ 



where a s — a(i] s ) and <f> s — <fi(r) s ). Both curvature and coupling are growing in this phase, which is expected to 
last until the time i] s when the curvature reaches the string scale and the background solution starts to deviate 
substantially from the lowest order solution. For recent ideas about how this may come about see [p|JTo|]. 

• The string phase lasts while ry s < r] < T)i. We assume that curvature stays high during the string phase. As in 
[plfl , we assume that the string phase ends when H (rji) / a(rji) ~ M s . We therefore implicitly assume that the 
string coupling at the end of the string phase is higher than it's value at the beginning of the string phase and 
that it does not decrease much during the string phase. But we do not assume any specific form of background 
evolution. Obviously, by doing that we give up on the possibility of obtaining any information about the 
produced particles during the string phase. We postpone such analysis until some better ideas and more reliable 
methods for handling the string phase are at hand (see, however, |^5| ). We do not lose the ability to compute 
the spectra produced during the dilaton-driven phase, as we explain later. We parametrized our ignorance about 
the string phase background, as in fl, by the ratios of the scale factor and the string coupling g(rf) = e^^^ 2 , 
at the beginning and end of the string phase z$ — a\/as and gi/gs, where g\ = e^ m " 2 and gs = e^ 8 " 2 . It 
turns out that also in the generalized setting we use, these two parameters are enough. We take the parameters 
to be in a range we consider reasonable, for example, z$ could be in the range 1 < zs < e 45 ~ 10 20 , to allow 
a large part of the observed Universe to originate in dilaton-driven phase, and gi/gs > 1- Note that unlike 
H or (lql we do not assume a specific background during the string phase. The time rji marks the start of 
the RD era as we explain below. Some other useful quantities that we will need are u>i, the frequency today, 
corresponding to the end of the string phase, estimated in [pj] , and the frequency u s = uJx/zs, the frequency 
today corresponding to the end of the dilaton-driven phase. 

Our definition of the string phase encompasses the whole duration in between the end of the dilaton-driven 
phase and the onset of RD phase, which may be quite complicated as in pQj. 
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The FRW RD phase is assumed to follow the string phase and last while r]i < r\ < n, 



— I = e^ 1 = const. (2.3) 



eq ■ 

a = ai ( — ] , — e^ 1 — const. (2.2) 

where n eq is the time of matter-radiation equality, a\ = 0(771) and (f>i — 4>(r}i) . Note that the dilaton is taken 
to be strictly constant, frozen at its value at n — r\\. 

• The FRW MD phase is assumed to last while n eq < n < r/o 

2 

where 770 is today's conformal time and a eq = a(n eq ). Note that the constant value of the dilaton is further 
assumed to be that of today. Some possible consequences of late time evolution of the dilaton were discussed in 

For more accurate results, we could have used a more sophisticated description of the RD and MD phases using 
entropy conservation, effective number of degrees of freedom and temperature, as was done in pi]] . 

We presented our model in the string frame, using string frame conformal time. Another favorite choice is the 
Einstein frame, in which the gravitational action takes the standard Einstein form. We continue to use the string 
frame exclusively throughout the paper. General arguments and many explicit examples have shown that physical 
quantities such as energy density or number density can be computed in any frame. 



B. The perturbation equation 

Deviations from homogeneity and isotropy are generated by quantum fluctuations around the homogeneous and 
isotropic background and then amplified by the accelerated expansion of the Universe. In addition to the dilaton 
and metric, string models contain many other fields that have trivial expectation values in our model of background 
evolution and do not affect the classical solutions, but they do fluctuate. We are interested in "low-mass" fields, which 
are either massless, or have masses much below the string scale, such as moduli, gauge bosons, and their superpartners. 

The action for each field's perturbation is obtained by expanding the 4 dimensional effective action of strings, 
which generically, for a tensor field of rank N, has the form i J d A x^f— <7e i *£ 2 (T Alli ... i/JJV ), where the parameter I in 
the dilaton prefactor is determined by the type of field. The second order covariant Lagrangian C 2 contains a kinetic 
term and possibly also a mass term. Setting the dilaton and metric at their background values, g^ v {n) = ^{rftn^, 
4>(j]) results in a quadratic action for each physical component of the perturbation, which we denote by ip, 

A=l J d V d 3 x a 2m e l * (^' 2 - (V^) 2 - M 2 a V) (2.4) 

where ' denotes d/dn, M is the mass of the perturbed field and m and I depend on the type of field. We ignore com- 
plication associated with gauge invariance and work directly with physical components. For example, the background 
metric has 1 co mponents of them 8 are non-physical, and each of the two physical components has an effective action 



of the form (2.4) 



As already mentioned, we are interested in particles for which the mass M is either vanishing, or much smaller 
than M s . In realistic models the mass may depend in a complicated way on the dilaton or other moduli. We will take 
M as an additional parameter and assume that it vanishes in the dilaton-driven phase and takes a constant value, 
much smaller than M s , from the start of the RD era and on. The parameters I and m take the following values in 
some particular cases. For the dilaton or graviton m = 1 and I = — 1, for the model independent axion m = I = 1, 
for perturbative gauge fields such as photons, m = and I = — 1, while for some non-perturbative gauge fields m = 
and I — 0, Ramond-Ramond axions and non-perturbative scalars have m = 1 and 1 = 0. In general, fields with more 
tensor indices tend to have smaller (more negative) m's. For example, the physical components of a 2-index tensor, 
whose field-strength is a 3-index tensor could have m = — 1. We do not know at the moment a general rule for allowed 
values of m and I, but it is obvious that many combinations are possible. 



The linearized equation of motion, satisfied by the field perturbation ip, derived from the action (2.4) is the following 



( A + M 2 ) + (l + (m - 1) (l - 1 /Vtj ) a~ V0' = (2.5) 
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where Atfj is the covariant Laplacian and we have used the fact that for our particul ar c ase = a(rj) 3+ ^. In 
Fourier space Aifi = a~ 2 Up'£ + 2Hip' k + fc 2- 0fe] and therefore the perturbation equation (2.5) is the following 

(2.6) 



i'l + \2H + (l + (m - 1) (l - 1/V3)) 0'] + {k 2 + M 2 a 2 ) ip k = 0. 



Using the transformation \ = S(r)) ip where S(rj) = a(rj) m g(r]) 1 and g(<fi) 



= P <t>/ 2 



, eq.(2.6) can be simplified, 



k 2 + M 2 a 2 - 



S 



Xk=0- 



(2.7) 



The function S(rj) was first introduced in |18|. 



C. Approximate forms and solutions of the perturbation equation 



Before trying to actually solve eq.(2.7) it is useful to look at approximate forms of the equation and their solutions 
in different physical situations. 

A few special moments exist in the life of a perturbation, first, 77 — > — 00 when it is essentially a free field perturbation 
in flat space-time. Then, the time rj ex when the rate of expansion H becomes too fast for a given perturbation and it 
"exits" the horizon and its amplitude freezes (see below). The perturbation remains "outside the horizon" until the 
"reentry" time rj re when the expansion rate H is again slow enough and the perturbation defreezes. We are interested 
in perturbations that exit during the dilaton-driven phase rj ex < i] s , and reenter during RD, rj re > 771 . We do assume 
that all perturbations stay outside the horizon during the string phase, using in an essential way that the curvature 
during the string phase stays high. 

In the range 77 < ?] e x, k 2 > S" / S cx H 2 , in the range -q ex < i] < r] re , k 2 + M 2 a 2 < S" / S cx H 2 and in the range 
i] > r/re, k 2 + M 2 a 2 > S"/S cx H 2 . Since we are interested in fields whose mass gets generated after the onset of RD 
and then have small masses, we observe that in the whole range —00 < r\ < ?/ re the mass term in the perturbation 
equation may be ignored. In effect, in this range of time, we may use the approximate form 

x£+(V-^)x fc = 0. (2.8) 

The general solution of the perturbation equation is a linear combination of two modes, one which is approximately 
constant outside the horizon and one which is generically time dependent outside the horizon. We understand the 
appearance of a constant mode as the freezing of the perturbation amplitude, since local physics is no longer active 
on such scales. The existence of the time dependent mode can be most easily understood in terms of a constant mode 
of the conjugate momentum of tp, II = S 2 {r])ip' . The amplitude of the conjugate momentum also freezes outside 
the horizon, since local physics is no longer active on such scales. This forces S 2 (r])ip' to be approximately constant 
leading to a "kinematical" time dependence of i[>. The perturbation for life is simplified if rewritten for \ n k = 5 _1 IIfe 

X^'fe + (k 2 { -^Pj X \ = (2.9) 

and it's solutions are Bessel functions, as for \k- 

In the RD Universe (MD Universe) a — airj/rji (a = a\r\ eq jr\\ (rj/ri elj ) 2 ), thus even if the S" / S is dominant at the 
beginning of the RD phase, then at later times M 2 a 2 + k 2 becomes dominant and therefore in the range 77 > ?/ re the 
perturbation equation may be approximated by 

X'fe' + (k 2 + M 2 a 2 ) X fe = 0. (2.10) 

The n umber density and energy density of the produced particles can be easily read off from the solutions of 
eq.( 2.10| ) in a standard way [|l6| 



We solve the perturbation equation using the following general method. We solve the early time equation with 
boundary conditions of normalized vacuum fluctuations. We assume that the constant mode of ip and the constant 
mode of II remain constant also while outside the horizon during the string phase, thus bridging the gap of unknown 
background evolution during the string phase. We then match the early time solutions to the late time solutions 
at 77 = rji when both the early time and late time equation are expected to be good approximations. For massive 
particles, an additional matching step is performed at 77 = 77 re . In all known cases our method reproduces results 
obtained by assuming specific background evolution during the string phase. 
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III. SOLUTIONS OF THE PERTURBATION EQUATION 



In this section we solve the perturbation equation, using the general procedure outlined in the previous section and 
obtain the necessary ingredients to assemble the spectra of produced particles. 



A. Early time solutions 

As explained previously, the approximate perturbation equation for early times is given by 



// 

Xk 



(3.1) 



In the range 77 < rj s , S(r}) — a ™ e l <t>sl 2 (rj/r/s) 1 ^ 2 ™ s where n s — + m) + 1 — m) /2, and in the range 77 > 771, 

S(jj) — a^e'^ 1 / 2 (77/771 ) 1 ^ 2 ~™ 1 where n x = i — m. Therefore, S" / S = (n 2 — i ) 77 2 , where 77 stands for either n s or m, 
and the perturbation equation (3.1) takes the following simple form 



Xl + (k 2 - (n 2 - i)77- 2 ) X k = 



(3.2) 



The solutions of this equation are Bessel functions. For the special values n s = ±i (corresponding for example 
to m = 0, I = 0), or 77,1 = ±2 (corresponding to m = 0,1), the equation simplifies and the solutions are simple 
exponentials. It is nevertheless possible to continue to use the general setup and substitute these special values at the 
end. 

Imposing that for krj — > —00 the solution corresponds to vacuum fluctuations ~^g e ~ lkr> ■> we obtain for 77 < r\ s 



while for for 77 > 771 we do not have any restrictions on the solution, 

Xk = Ay/rjJ n (kr)) + B^/r]Y n (kr)). 
For 77 < 77s, for perturbations outside the horizon, |A;?7| < 1, Xk can be approximated by 



Xk 



/7F (1 — i cot (n s 7r)) 
2 r(l + n«) 



1 - 



(fcT?) 2 



2(2 + 2n s ) 



. ^/tt esc (ji s tt) 



1 - 



(fcT,) 2 



2(2 - 2n s 



while for 77 > 771 and krj < 1, Xk can be approximated by 



(1 - 7 cot (n l7 r)) 

r (1 + 771) 



1 - 



(kr,) 2 



2(2 + 2ni) 



csc(mTr) 

F(l - 77l) 



Therefore, in the range rj < rj s the solution is approximately given by 

2-i-n. 



-r(-77 s ) as T m e-^/ 2 (fc s fc)"= 77 2 ™= 



nn,-l / 7, 

+ ^=r(77>r n e-^ /2/ " 

V7TK. 



0^) 2 



1 - 



(fc??) 2 
2(2-2t7i) 



2(2 + 2n s ) 

(fcT,) 2 



2(2 - 2n s 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



where we have used the relations r(n)r(l — n) = n/ sin(77-7r) and r(— n)T(l + n) = 7r|l — icot(7rn)|, and defined 
k s = l/|?7s|. In the range 77 > 771 we obtain 



out a 

k — A 



1 



27rVfci 



T (-m) ar m e- i0l/2 (fcifc)" 1 77 2 " 1 



{krif 



D- 



^2™ l r(77 1 )ar m e^ l/2 (£■ 



TX\[k\ 



2(2 + 2m) 
(fc?7) 2 
2(2 -2ni) 



(3.8) 



G 



where we have used the relations T(n)r(l — n) — ir/ sin(n7r) and r(— n)F(l + n) = tt|1 — icot(7rn)| and defined 
ki = 1/ryi. Note that we have kept, for later use, also the next to leading terms. The corrections are of order (kr]) 2 , 
which is a very small quantity, showing that the approximation we use is quite good. 
The perturbation equation for \™ 1S the following 



Xl"+\k 



(s- 1 )" 
s- 1 



xl = 0, 



(3.9) 



and it's solutions are Bessel functions, as for Xk- We normalize the solution at early times to vacuum fluctuations, 
_/| e -ifc?7 ( no t e the \fk instead of \ j\fk in the normalization of the momentum). Substituting the explicit expressions 

for S(r]) = a m {r))e 1 ^ 2 , where a(rj) and <p(r)) are given in eq.(^J), eq.( |3.9[ ) takes a form similar to eq. (|3.2| ), except 
that instead of n s appears and instead of Hi appears nj, where = n s - 1 = (V3(Z + m) - 1 - to) /2 and 
n\ — ni — 1 = — | — to. The solution of eq.( |3.9[ ) for r\ < r] s is the following 



(3.10) 



and for r] > r]i 



Xl = A^}J n *(kr]) + By/fjY n «(kTt). 



(3.11) 



For the special values = ±5 and n\ = ±^ the solution reduces to a sum of exponentials. It is nevertheless possible 
to use the general formulae for these special values as well. 

For 77 < r] s , for perturbations outside the horizon, \kr)\ < 1, x n can be approximated by 



w r(-nj). _ 



(fcT?) 2 



2(2 + 2nJ) 



(fcr,) 2 



2(2 - 2<) 



(3.12) 



while for rj > r/i, kr) ^ 1, x n can be approximated by 



(kr)) 2 



2(2 + 2n\) 

Since for r] < r) s , II = S'x 7 ' = afe 1 **' 2 (r)/r} s ) X,2 ~ ns 



£K) /fcr? 



(fc??) 2 



2(2 - 2<) 



(3.13) 



IK 



"m^" — 1 



-r (-nj) oj'e'*'/ 2 (fc s ) 1/2 



(fcry) 2 



2(2 + 2toJ) 



-r «) a^ /2 fc (fc s r 1/2 (fc s fc)- a+1 rr 2r ' 



(fcr,) 2 



2(2- 2nl 



(3.14) 



and since for 77 > 771, II = 5x 7r = a^e'* 1 / 2 (77/771 

2 -< 



,1/a-m ^ 



/ 1. \ ni 

IT',"" = = J 4 7r r(-n^)a r 1 "e^ 1 / 2 fc- 1 fc 1 1/2 f- 

i— B„T «) are^ l/2 fc- 1/2 (fcxfc)""^ 1 77- 2r 



(kr]) 2 



2(2 + 2<)_ 
(kr]) 



1 - 



2(2 -2raj) 



(3.15) 



ing t/ie gap: The constant mode of ip 



We will compute the contribution to the out perturbation from the constant mode of ip m . From the second term 
of eq.(3/7) we obtain the the constant mode of ip m 



4i 



nFk 



-_2 n °- x Y (n s ) a7 ro e-^»/ 2 



(3.16) 
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For 77 > 771 , ip out is given in eq.(B.8). We assume that the constant mode of 4> is indeed constant during the string 
phase (and, in general, outside the horizon), as explained in the previous section, and use this to connect the two 
solutions by matching them and their first derivatives at r) = 771. This matching procedure is expected to be quite 
accurate because we match the solutions at a point in time at which both solutions are well outside the horizon. We 
pres ent the m atch ing equations once explicitly, to emphasize the importance of keeping the next to leading terms in 
eqs.Q and ©. 



1 



2 ni 



B 



-T(n 1 )a^ m e- 14 ' 1 / 2 



B- 



(km) 2 



rji 2(2-2711)71-^ 

o-ni 

A- (-m) a^e' 1 ^ 1 ' 2 



A 

k s 

hj 

r (m) ai 



A- 



2 -m 



m -ltp 1 /2 



T(~ ni )a^ n e- 1 ^/ 2 (hk) ni r? 2 ™ 1 , 

k 
k[ 



n 1 \ n l 2ni 

(kik) r\ x 



The solution of eq.( 3.17 ) is 
A = 



2 n s + ni -2^ r(yts) +m+ni ^ 

m(l-m) T(- ni ) Zs 



B 



r\n s —Til — 1 



-£ (n s ) 1/2+m-m / gs_ 



91 

-1 



k 

ni—n 



2—ni — n s 



o 0i/2 



gs 



(3.17) 

(3.18) 
(3.19) 

e^/ 2 . Note that had we not 



To obtain eqs.(3.1£ ),(3.19) we have used the definitions zs = ai/as, gi 
kept the next to leading terms in eq.(3/7) we would have obtain ed A = 0. 

To summarize, plugging in the resulting coefficients from eqs. fl3T8l ) and fl3T9| ) into eq.(§J) we obtain 



' 2— r ( „, )( .r"' 9 ;'4™(i 



1 



2™ s r(n s ) m , _j 
-u y s z s 



'irks ni(l — m) 



2-n s 



ftjo 



(3.20) 



The out solution ( 3.20| ) may be put, using the the relation /cs 2 s = fci, in the form ip™* = Ci 



i + C 2 & (M) 



\2ni 



where C2 is a numerical coefficient of order 1. At 77 = rji, the ratio of the two terms is approximately (krji) 2 <C 1, 
because kir/i — 1. However, at later times, using eq.([0|), the general solution is given by ip% ut ~ ^4 S~ 1 {rj)^/rjJ ni (kr}) + 
B S~ 1 (ri)y/rjY n *(kr)). A good time to look at is near 77 re . At that time the magnitude of both Bessel functions is 
generically the same, and to decide on the leading contribution we look at the ratio of the coefficients A/B ~ 

, 2-2ni 

Therefore the leading contribution at much later times depends on ri\. For ni > 1, A 3> B and the 



leading contribution will be from the time depen dent term in eq.( 3.2C ), and for ri\ < 1, B 3> A and the leading 
contribution will be from the constant term in eq.( |3.20 ). In conclusion, for later times, 



1>%*~AS-\r,)y/rjJnM 



ni > 1 
m < 1, 



(3.21) 
(3.22) 



where A and B are given by eq.(fT8|) and eq.([yj| 



2. Bridging the gap: The time dependent mode of if) 
We will compute the contribution to the out perturbation from the time dependent mode of ip m . The time dependent 



mode of ij) can be extracted from the first term in eq.(3.7) 

1 



2V7rfc7 



Since for 77 < r/ s , S 2 (r)) = a 2m e 1 ^ (77/rfc) 



l-2n s 



T{-n s )a- m e~ 1 ^' 2 (fc s fc)" s V 2 



, we obtain 



(3.23) 
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W s 2 



T(-n 6 



la m e /0 s /2 



const. 



(3.24) 



We can see explicitly that the time dependent mode of ip corresponds to the constant mode of the canonically conjugate 
momentum II. 

We compute the contribution from the time dependent mode of ip using the constant mode of II, assuming that it is 
constant throughout the string phase (and, in general, outside the horizon), as explained in the previous section, and 
use this to connect the two solutions by matching them and their first derivatives at r\ = rji . This matching procedure 
is expected to be quite accurate because we match the solutions at a point in time at which both solutions are well 
outside the horizon. 



The solutions are given in eq.(B.14) and (3.15) and their matching yields 



An 



£>7T 



v ^ 2 m-n.-ir(l-n.) -i/2-r.+m f 9s V fe ( k_ 



i3ni 



r(i-m 

~ 2 r(l-n 



,9i 



-ni+n s 



— 1/2— m— Tii 



m(i - 2n x ) r(m - i) 



ni -\-n s 



where we have used the definitions of and n\. Substituting into eq.(3.15) we obtain 



nr* 



2 i-«i 



A v T(l-n 1 )aTe l ^ 2 k- x k 1 1 /2 (^- 

V fci 



-B„T (m - 1) a\ n e l ^l 2 k- 1/2 (k 1 ky ni+1 V 2 - 2ni 



(3.25) 
(3.26) 

(3.27) 



with An and B^ given in eq.( ft.25| ) and ( |3.26 ). 

As for the case of the constant mode of ip, we may compare the relative strength of the two terms. At n = r/i , the 
relative strength of the two terms is approximately (krji) 2 -C 1, because kirji — 1. However, at later times, using 
eq.(3.11), the general solution is given by U^ ut ~ A^ S{n)^/r]J n ^{kvi) + B n S \n) y/rjY n * (krj) . A good time to look at 
is near r\ re . At that time the magnitude of both Bessel functions is generically the same and to decide on the leading 

/ , \ -2th 

contribution we look at the ratio of the coefficients A^/B^ ~ ( J . Therefore the leading contribution at much 

later t imes depends on n\. For ri\ > 0, A n 3> B^ and the leading contribution will be from the time depen dent term 
in eq.( 3.27 ), and for n\ < 0, B n 3> A T and the leading contribution will be from the constant term in eq.( 3.27 ). In 
conclusion, for later times, 



* S(n) A^J nl {krf) for m > 

* S(n) B„^rjY n -«(kr)) for m < 

where A^ and B^ are given in eqs.(3.25), (|3.26) . 



(3.28) 



B. Late time solutions and spectrum: The massless approximation 

The perturbation equation in the range r\ > r\ re is x'l + f^ 2 + M 2 a 2 — ^-j \k — 0. Because we cannot solve this 
equation in a closed form for the general case, we use an approximate form of the equation to obtain approximate 
solutions. We distinguish between two cases, k > Ma re - the "massless approximation"; k < Ma re - the "massive 
approximation", where a re — a{rj re ). 

In the massless approximation the perturbation equation is x'k + f^ 2 ~ Xfe = 0> which is the same equation we 

solved already, so there's no need for any additional steps. All we need to do is to take the late time solutions and 
interpret them. In a standard way, we compute from the solution the number of produced particles at a given wave 
number k, iVfe, and the ratio of the energy density per octave of the produced particles to the critical energy density 
-.t^- . For massless particles <Kl/ {dhiuj) = w 4 N lu /(M p Hq) 2 , where Hq is today's Hubble parameter, M p is the Planck 
mass and u> — k/a(n) p4| . 

We can further simplify the denominator (M p Hq) 2 . For that we switch momentarily to using cosmic time. We may 
reexpress the relation Hq = H eq ^-, as H Q = Hij^^-, where H 1 is the Hubble parameter at i l5 the onset of RD. In MD 

2/3 , s-3/2 _ , N l/2 



a(t) = a eq and therefore ^ = (^j . In RD a eq = a x (j^j , then ^ = C^A . Combining these 

2 3/2 2l/2 

relations, we obtain H Q = Hi (^A (f^ . Since we assumed that Hi = M s then H = M s (f^-) 
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Now, M p = M s g 1 , M s — (Jiao/ai and u>% — lo s zs, and therefore M s ^ = lo s zs- Combining these relations we arrive 
at our final destination 



w^I^V 1 . (3-29) 



where z eq = a^ja e 



Substituting eq.(3.29) into the expression for the energy density we obtain the following expression 

eq 



dn/(d\ncj) = z-^g'i I — ) N u . (3.30) 



The factor z^ 1 , the redshift from the time of matter-radiation equality and today, appears because the energy density 
of massless particles redshits faster than that of matter. 

We may also discuss spectra of massive particles in the massless approximation. Although this may sound a little 
strange, such an approximation does make sense, because it is made at rj re , when the effects of the mass can be 
negligible. The only difference between massive particles and massless particles in the massless approximation is 
in their energy density. The solutions near r\ re are approximately the same, and therefore the number of produced 
particles is the same as in the massless case. However, at much later times, the contribution from the mass may 
dominate the energy density. We can then compute the ratio of the energy density per octave of the produced 
particles to the critical energy density ^jf^j- For massive particles dfl/(d\ncj) = VA/ 2 + lo 2 w 3 N u / ( RI v Hq) 2 , where 



Hq is today's Hubble parameter, M p is the Planck mass and u> = k/a(rf). Using cq.( 3.2E ) we may show that 



dSl/(d\nu) = z'J-z^gl v ^ — I — ) N u . (3.31) 



, \JM 2 + a- / X 



An important difference compared with the massless spectrum is that one power of uj/lo s is replaced by ^ M ^+ u)2 
(for u) < M). Because the energy density of massive particles redshits slower than that of radiation, two enhancement 
factors will appear compared with the case of massless particles. One enhancement factor accounts for the period 
between reentry and matter-radiation equality, in which the relative energy density of massive particles grows, while 
that of massless particles remains constant. Another enhancement factor appears because the relative energy density 
of massive particles stays constant during the matter dominated era, while that of massless particles decreases. It is 
harder to see this factor in the equations because z eq also appears in the expressions for massive particles. However, 
it's appearance is artificial. The combination M/(z eq u s ) = M / 'to s (rj eq ) actually does not depend on the redshift from 
matter-radiation equality epoch until today. 



1. Massless spectrum: contribution of the constant mode of ip 



The approximate late time solutions of the perturbation equation are given in eqs.( 3.21 ), ( 3.22 ). Using the asymp- 
totic expansion of Besscl functions, the solutions (either during RD or MD) approach, for krj > 1 



where 



and 



iPk = S- 1 ^=(a k e- ikr >+(3 k e ik *>) 
V 2fc 



Kl = Wk\ 



nn s — ni 



r(n.) 



Kl = \Pk\ 



)n s +ni — 1 



m(i - ni) 



r( ni ) 
r(« 



r(-m) 



l/2+m-ni ( 9S_ 
91 



-3/2+m+m ( 9S_ 

9i 



for n\ < 1 



-/ 



2— ni — n B 



for n\ > 1 



(3.32) 



(3.33) 



(3.34) 



The number of massless particles Nk, can be expressed in the standard form Nk = \Pk\ ■ Using eq.(3.30), ri\ 
and using d£l^ to denote the contribution from the constant mode of ip w e arrive at the final expressions 
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da 



din uj 



i)2n s +2m 



dint 



ni<l 



ni>l 



r 2 K) 



b-2m-2n s 



7+2m-2n s 



(3.35) 
(3.36) 



2. Massless spectrum: contribution of the time dependent mode of ip 



As usual, we use the constant mode of II to evaluate the contribution of the time dependent mode of ip. The 
approximate late time solutions of the perturbation equation are given in eq.( |3.28 ). Using the asymptotic expansion 
of Bessel functions, the solutions (either during RD or MD) approach, for kr\ > 1 



11* 



S-j= (a k e 



ikrj 



ikrj\ 



where 



and 



T{-nl) 



-l/2-m+ni ( 9S_ 
91 



-ni+n s 



for m > 



2"3ni — n s — 1 

|ni(l-2m)| 



r(-nj) 



r(nf) 



-1/2-m-m f 3s 
3i 



for ni < 0. 



(3.37) 



(3.38) 



(3.39) 



The number of massless particles Nk of a given wavenumber k, can be expressed as Nk = |j9fc| ■ Note that we are 
calculating the number of particles directly from the solution for Ilk- We could have also used the relation II = S 2 ip', 
use the stan dard definition of Nk from the coefficients in the solution for ip and obtain the same result. 

Using eq.(3.3C), n\ = \ — m, and using dil^ to denote the contribution from the time dependent mode of ip (the 
constant mode of II) we arrive at the final expressions 



dQ T , 
din lj 

dim 



T^2 /i \ / \ 2/ / \ 3+2m+2n B 

_ n l-2m-2n e T U ~ n sj -1 2 -im-i 9S\ ( V 

I*(m + §)*-»** fl UJ U. 

-6m-2ri s -5 r 2 /i „ \ /„ \ 2 ' / , , \ 5-2m+2n s 
1 (l-n s ) _j 2 -6/3SI / W 

,2^ _ 1^2 r2 I _ ^ 2 eg 01 Z S 



rii<0 



m 2 (m- i) 2 T 2 (-| 



TO 



.91 



(3.40) 
(3.41) 



3. Massive spectrum: contribution of the constant mode of ip 

The final result for the contribution of the constant mode of ip to the spectrum of massive particles (recall that we 
are using the massless approximation, defined at the beginning of this section) maybe read off from the final result for 
massless particles, eqs.( 3.35 ) and ( ft.36| ), by making the substitution of one power of uj/uj s by one power of xME+i^ 



da 



dlnu 



= 2 Z 



,+2m-l ■* 2 ( n s) ,-l„2.4m-4 ( 9S 



da 



dlnuj 



ni<l 



ni>l 



)2n s —2m— 1 



r 2 (±-m) 
r 2 (n.) 



z e 9 3i z s — — 



.9i 



4-2m-2n 3 



Vm 2 



.-X-2_-6 f 3S 



(i-m 2 ) 2 T 2 (m-i) 



Z e?5l% [ — } \ — 



-21 



e+2m-2n s 



r2 



(3.42) 
(3.43) 
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4- Massive spectrum: contribution of the time dependent mode of ip 



The final result for the contribution of the constant mode of ip to the spectrum of massive particles (recall that we 
are using the massless approximation, defined at the beginning of this section) maybe read off from the final result for 
massless particles, eqs.( 3.4C ) and ( p\4l| ), by making the substitution of one power of uj/uj s by one power of ^ M '^+ ujZ ; 



din a; 



= 2 



l-2m-2n. 



r 2 (1 - n s ) 



dhxoj 



rn>0 



ni<0 



r 2 (m + 

2 -6m-2n,-5 r 2 (1 - n s ) 
m 2 (m — i) 2 r 2 (— i — m 



-1 2 -4m-4 
9l z S 



< <l 



-eg 



21 



9i z s — — 



2+2m+2ra s 



4-2m+2n s 



VM 2 



Vm 2 



(3.44) 
(3.45) 



C. Late time solutions and spectrum: The massive approximation 



If the produced particles are massive, and at f] re , Ma re > k, we are unable to solve the late time equation exactly, 
and therefore need to resort to another matching procedure. We choose to perform the matching at r\ = rj re , where 
k 2 + M 2 a 2 e ~ H 2 . We match a solution from outside the horizon to a massive solution inside. This is not an ideal 
match, but both are expected to be fairly accurate at rj re . Note that the solutions outside the horizon were obtained 
using a "massless" equation, however outside the horizon the dominant term is S" / S anyhow, and therefore no large 
inaccuracies should be induced. We then need to know whether r] re occurs during RD or MD. If Ma eq > H(r} eq ) 
then the the reentry time of all perturbation will be at earlier times. Therefore, for M > 10~ 27 eV, the reentry time 
is during RD. Since we do not expect any particles to have masses smaller than 10~ 27 eV, we continue under the 
assumption that r) re < r\ eq . We have performed the calculations for the case r\ re > r) eq , however, we will not present 
them in this paper. 

In a standard way, we compute from the solution the number of produced particles at a given wave number k, N^. 
We then compute the ratio of the energy density per octave of the produced particles to the critical energy density 
''" As already shown, for massive particles dfi/(dlnw) is given by eq.(3.31). 



d In i 



The equation of motion in the massive approximation is 



X £ + M 2 a 2 Xk = 



(3.46) 



and it's solutions in the WKB approximation are given by 



xt 



I cxp (±i [ Ma(n')dn' 

Ia(n) y \ J 



(3.47) 



where for the RD epoch (77 < r/ eq ) 



xt 



1 / 1 „ 

= exp ±i— Maikxn 



(3.48) 



Since ip = S 1 \ the general solution for the massive perturbation for RD is the following, 



massive 



1 



y /2Maik 1 rj 



ci cxp I -i-Maikirj 2 j + (3 k exp (i-Maikirf 



(3.49) 



We will also need the general solution for the conjugate momentum, which can either be derived form the equation 
for x* or from the relation II = S 2 ip' , both yielding the approximate solution 



n massive 
k 



y/2Maik 1 r)S(r)) 



n f exp J -i^Ma 1 k 1 i] 2 j + (31 exp (i^Maikirj 2 



(3.50) 



Instead of using only ip at late times, we will again use II, making sure that we do not double count the number of 
produced particles. 
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1. Massive spectrum: The contribution of the constant mode of tp 



A s in t he mas sless case, we assume that H 2 > k 2 + M 2 a 2 , for all k, therefore solutions for r\ > r\\ are given in 
eqs.( 3.21 ), ( [3.22 ). We reproduce here their approximate forms, 



m(l-ni) r(-ni) 



fjl 



-I / i \ 2— m —rii 



y/rjJ ni (k-q) for > 1 



,-m-l£K) y 2n 

YW) Zs 



fJS 



-/ 



~ ) I 7T ) VvJ-m {hrf) for m < 1 



(3.51) 



(3.52) 



In the case of the massive approximation these solutions are correct only for H 2 > k 2 + M 2 a 2 . Therefore for kr\ < 1 
(and M < k/a) the leading terms are 

2-n 3 



1 2^- 2 r(n s )_ m ^ 2(ni _ 1) 



/7rfc s ni(l - rix) 



for the time dependent mode, and 

V>m<i 



3s 



(3.53) 



(3.54) 



for the constant mode. 
(a) The case n\ > 1. 

In the case m > 1, the ap prox imate solution for r\ < rj re is given in cq.( 3.53| ), and for 77 > r\ Te we obtained the 
approximate solution in eq.( |3.49| ). We connect the solutions at r\ = rj re by matching them and their first derivatives, 
V^'O/Ve) = ^ assive (r] re ) and tp^iVre)' = ip^ assive (rjre)' ■ Using the relations H(r] re ) = Ma(r? re ) and t?" 1 ee fc re = 

1 /2 

(Mai^i) 7 , the results of the matching give 



|3m + z| 2 



n„-3 



l r ("«)l „-i ( 9s 

H |m(i-m)| s 



91 



-/ 



! " . s . an 

|3m-M| 2"°~ 3 |r(n 8 )| 
0r |m(l - ni)| 



Since k re /k s = (Mai/ki) zg, and since we assume as in plfl fci/ai = M s , we obtain 

-2 / \ 2— 71 



,TO-3/2 f 3£ 
31 



M 

Mi 



ni/2 



Using eq.(3.31), and the definition of n±, we obtain the final result 



d\nu) 



|3m + i| 2 2 2n °- 6 



ni>l 



(J-m>): 



-r 2 K)z-; 5 2 4™- 7 (^ 

\3i 



VM 2 + w 2 / u> 



7-2n s 



M 

Ml 



m-1/2 



(3.55) 



(3.56) 



(3.57) 



(b) The case n± < 1. 

In the case ni < 1 we obtained the approximate form of the solution for rj < rj re in cq.( |3.54 ). and for r\ > rj re 
we obtained the approximate form of the solution in eq.( |3.49D . We connect the solutions at 77 = rj re by matching 
them and their first derivatives, ip^iVre) = yb]? assive {rp^j&nd ^ ut (T] re )' = ^ asslve {r lre ) 1 . Using the relations 



H(r) re ) = Ma(r] re ) and r/ rf } = k re = (Afaifci) 1 ^ 2 , the results of the matching give 



Since k re /k s = (Mai/ki) 1 ^ 2 zg and since we assume as in |2l[| k\/a\ — M s we obtain 



+ m + i\ 



1/3* 



+ m + i 



l r K)k' 



1+1/2 f ffS 

3i 



Af 



ni/2 



Using eq.(3.31), and the definition of ni, we obtain the final result 



dlnuj 



+ m + i\ 



ni<l 



; 2 2 - 4 r 2 (n.) z- 5 2 4 m - 3 (f ^ VM2+ ^ ( " 



M 

Ml 



1/2-r 



(3.58) 



(3.59) 



(3.60) 



Note that the spectrum is actually partially amplified because the relative density of massive particles during RD is 
growing as the scale factor. 
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2. Massive spectrum: The contribution of the time dependent mode of ip 



As in the massless approximation, we assume that Hf > k 2 + M 2 a 2 , for all k. We use our method of calculating the 
contribution of the time dependent mode of ip by using it's correspondence to the constant mode of II. The solutions 
are given in cq.( |3.28| ). We reproduce here their approximate forms, 



U° k ut = S(rj)^ ni ~ 



-I r (-Q -2m ( 9S_ 

\9i 



-ni+n s 



n x (l-2ni) T(n1) 



l / \ ni+n 



k %Jr\J n \ (krf) for n\ > 



(3.61) 
(3.62) 



In the case of massive particles these solutions are correct only for H 2 > k 2 + M 2 a 2 . Therefore for kr\ < 1 (and 
M < k/a) the leading terms are 



n? 



1-nl 



r(-<)« 



for n\ > 0, 



and 



o4ni — n s — 3 

nr* = \/^kl^ —Y{-nl)a n s l g l s z s 1+2m 



' ni(l — 2ni) 



l+n s — 2m 



(kf]) 2m+1 for m < 0. 



(3.63) 



(3.64) 



(a) The case ri\ > 0. 

In the case ri\ > 0, the approximate solution for rj < r/ re is given in eq.(3.63) and for r\ > r\ re it is given in eq.( |3.50 ). 
We connect the solutions at r\ = r\ re by matching them and their first derivatives, W° k ut (r\ re ) = Ji™ assme (jj re j and 
n° k ut (r] re y = tt™ assive (r] re y, and obtain 



|/3fc| = \\+m + i\2- 2 -<^\T{-n*)\ 



-2m ( 9S_ 
.91 



Using k re /k s — (Mai/ki) 1 ' 2 z$ and since we assume as in |21j k\ja\ — M s 



1/3*1 = |±+m + i|2- 



^\T(-n:)\z- 



m-1/2 ( 9S_ 

9i 



M 
Ml 



(3.65) 



-ni/2 



Using eq.( 3.31 ), and the definitions of and n-y, we get the final result 

2; 



dine 



ni>0 



5s V VM 2 + uo 2 ( u 
91 



3+2n s 



M 

Ml 



m-l/2 



(3.66) 



(3.67) 



(b) The case ri\ < 0. 

In the case n\ < 0, the approximate solution for 77 < r\ re is given in eq.( 3.64 ) and for 77 > r\ Te it is given in eq.( |3.50 ). 
We connect the solutions at rj = r\ re by matching them and their first derivatives, U^. ut (r] re ) — U k nassvve (r] re ) and 
Uf, ut (r] re y — U k nass ' lve (ri re y . Using k re /k s — {Mai/k s ) X ^ 2 z s , and since we assume as in plj k±/ai — M s , we obtain 



\ph\ = \-\-m + i\ 



4ni — n s —4 



ni(l — 2n\) 



|r(-nj)|zs 



5/2-m f ffs 
3i 



2+n s 



M 

M 



ni/2-1 



(3.68) 



Using eq.(3.31), and the definitions of rig and ni, we obtain the final result 



c?ln a; 



-6-2n s -8m 



— m + l\ 7T 



ni<0 



m 2 (i — to) 2 



-1 2 -9-2m / 9S_ 
z eq 9l z S I 



r 2 (l -ra s ) X 

7+2ra a 



2 + bJ 2 ( U) 



M 
Ml 



-m-3/2 



(3.69) 
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IV. SPECTRA OF PRODUCED PARTICLES 



In this section we collect the different contributions that we calculated in the previous section and put them together 
into complete spectra for massless and massive particles. We use this opportunity to recall the parameters that appear 
in the expressions for the spectra. The spectra we obtain depend on two parameters, the total redshift during the 
string phase z$ — ai/ag, and the ratio of string coupling at the beginning and end of string phase gs/gi- The 
coupling gi is considered a known number, the value of the string coupling today, and so is the string mass M s . The 
frequency oj s is expressed as oj s = W\jzs where u>\ is taken to be known |21|. Also appearing are the redshift since 
matter-radiation equality z eq . We recall also the definition, n s = U/3(l + m) + 1 — ml /2. All spectra are valid for 
lu < lo Si unless otherwise stated. We present the resulting spectra in their original form, but also in a more symmetric 
form in the spirit of |p.8| , to help in the possible uncovering and understanding of some underlying symmetry. 



A. The spectrum of massless particles 
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+ [jjtzj] ni> Q- The relevant expressions appear in eq.( 3.36 ) and ( 3.40 ) 
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The leading contribution is determined by the ratio of the two terms which is approximately z^r 4 "^^-) 4 ' , if ^ < 1 
and ^5 > 1 then the second term will dominate. The frequency dependence of the second (and leading) term in 
cq.(4.1) is determined by the index n = 3 + 2m + 2n s = 4 + m + V3(l + m). So, for example, for a perturbative 
dilaton dependence / = —1, n = 4 — V3 + (1 + ^/3)m. Therefore, all spectra with I = —1 and m < — 1 will be 
decreasing spectra. In particular, the case I = — 1 and m = — 1 corresponding to the usual antisymmetric tensor will 
have n = 3 — 2\/3 ~ —0.46. More negative values of m will give more sharply decreasing spectra. 

We may present eq.(i.l) in a more symmetric form suggested in p8[. We do this to allow for a future study to 



reveal underlying symmetry. 
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where Af/^ and AY m can be read off eq.(4.1) 
(b) The case 1 > ni > (\ > m > -\). 
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. The relevant expressions are given in eq. (^35| ) and ( pp| ) 
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The leading contribution is determined by the ratio of the two terms which is approximately z s 8m (^-) 41 , if ^ < 1 
and zs 1 then the dominant term is determined by the sign of m. The case m = 0, I = — 1 corresponds to photons, 
studied in |l2|. In that case, the second term dominates because of the factor (gi/gs) 4 - It is the only term computed 
in [l2] ] . The index determining the frequency dependence is n = 4 — \/3. It is comforting to observe that our method 
of calculation, not relying on specific background evolution during the string phase gives the same results (where 
available) as a calculation relying on some specific background evolution. 
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We may represent eq.(4.3) in a more symmetric form. We do this to allow for a future study to reveal underlying 
symmetry. 
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where N( m and Af* m can be read off eq. ( |4.3| ) 

(c) The case n\ < (m > \). 



In this case 
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[ ,, a 1 . The relevant expressions are given in eq. (f]35j) and ( ^4l| ) 
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(4.5) 



The leading contribution is determined by the ratio of the two terms which is approximately 



y — Am— 2(gs \4l 



/SS)4( if 



^ < 1 and zj > 1 then the first term will dominate. The frequency dependence of the first (and leading) term in 



eq.(4.5) is determined by the index n — 5 — 2m — 2n s — 4 — m — V3(l + m). So, for example, for a perturbative 
dilaton dependence / = — 1, n = 4+ vo - (1 + v3)m. Therefore, all spectra with I = — 1 and m < 1 will be increasing 
spectra. In particular, the case i = — 1 and m = 1 corresponding to gravitons and dilatons will have n — 3 as in [[L7p|] 
(we have used an approximation in which the logarithmic factors are missing). Again, we observe that our method 
of calculation, not relying on specific background evolution during the string phase gives the same results (where 
available) as a calculation relying on some specific background evolution. 

More negative values of m will correspondingly give more sharply increasing spectra. The spectrum of the model 
independent axion, for which I = 1, m = 1, will have an index n — 4 — 2\/3 ~ —0.46. Note that this is exactly the 
index of the leading contribution to the spectrum of the antisymmetric tensor (a fact known to the authors of ]l9|). 
In this class of generically increasing spectra the axion is standing out with it's decreasing spectrum. 

We may represent eq.(4.5) in a more symmetric form. We do this to allow for a future study to reveal underlying 
symmetry. 
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where N?U and A\ *1 can be read off eq. (|4.5|) 



B. The spectrum of massive particles 

1. The discontinuity in the spectrum of the massive particles 

The spectrum we computed in section (3.C) was that of massive particles that enter the horizon as massive particles, 
so at r] re , k/a{rj re ) < M. The spectrum of massive particles computed in section (3.B), is valid for massive particles 
if at rj rei k/a{rj re ) > M. In the approximation scheme that we are using, there will be a discontinuity in the slope 
of the spectrum of massive particles, at a wavenumber k m , where k m is defined by the condition k m /a{rj m ) = M. 
The moment r\ m is the (conformal) time where the approximation k/a(r] m ) > M at reentry is no longer correct. To 
compute k m we switch momentarily to cosmic time. First, since t m is during RD (see below), it can be computed 
from M = H(t m ) = Hiti/t m . Since a(t m ) depends on whether t m is bigger or smaller then t eq , we get that if 
tm < teq (he. M > H(t eq ) — lCP^ev), all the particles will get into the horizon in the RD period. Therefore 
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(t m ) = ai(t m /ti) 1 / 2 = a^Hi/M) 1 / 2 and since H% = kx/ai we obtain a{t rn ) = (fciai/M) 1 / 2 , 
Secause ^ = V and fa/ai = M s . V = (M ° 1 f fcl)1/2 . Therefore 
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The 
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discontinuity of the slope will be at a frequency uj m = lo\ I -p- J . At that frequency 
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If w s > Af , there will be an additional change in the slope at uo = M, abov 
of massive particles and massless particles become approximately the same. 



(4.8) 
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2. Spectrum of massive particles 
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For uj m < uj < lo s the relevant expressions are given in eq.(3.43) and (3.44). 
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For a; < w m , the relevant expressions are given in eq.(3.57) and (3.67) 
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(4.10) 



that the spectrum is continuous, up to numerical factors of order 1. The 



At uj — oj m we may check, using eq. 
spectral index, however, jumps at oj m . 

I f uj k > M, then for the range M < oj < oj s the spectrum is approximately the same as the massless spectrum, 



eq.(4.1 



As in the massless case we may determine the leading contribution to the spectrum. It is determined by the ratio 
of the terms which is approximately z 2 s ~ Am , if ^ < 1 and z$ ^> 1 then the second term will dominate. For 

example, for I = —1 and m = —1 the spectral index will be n = 2 — 2^/3 ~ —1.46 for w > uj m and n = 5 — 2^3 ~ +1.53 
for w < uj m . I f oj s > M, then for M < lu < uj s the index is n = 3 - 2\/3 ~ -0.46. 



Eqs. ( |4-S| ) , (4.1C) may be presented in a more symmetric, perhaps related to some symmetry of the underlying 
physics. We present here only the symmetric form of eq.(4.1C). 



dn 

dlnuj 



Ml 



1 2 -6 
l,m z eq 9l z S 



VM 2 + OJ 2 



-21 

9S \ _2m-l 



X < Km \ — } s 











2-2n s 




) 





21 



-2m+l 



-2+2n s 



(4.11) 



where M.\ m and B\ m are numerical coefficients of order one which could be read off eq.( 4.10 ) 
(b) The case 1 > m > (— | < m < |) . 
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For u> m < uj < uj s the relevant expressions are given in eq.(3.42) and (3.44). 
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For u> < u> rn , the relevant expressions are given in eq.(3.60) and (3.67). 
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At uj — uj m we may check, using eq.(l.S) that the spectrum is continuous, up to numerical factors of order 1. 

If uj s > M, then for the range M < uj < uj s the spectrum is approximately the same as the massless spectrum, 
eq.Q. 

The determination of the leading contribution is more complicated now, because the ratio of the two terms in each 
range is different and involves also the ratio j^-. We postpone such analysis to a more detailed study of specific cases. 

Eqs.( [t~12" ), ( 4.13 ) may be presented in a more symm etric, perhaps related to some symmetry of the underlying 
physics. We present here only the symmetric form of eq.( 4.13| ). 
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where M 1 ^ and Bf^ are numerical coefficients of order one which could be read off eq.( 4.13 ). 
(c) The case ri\ < (m > i) . 



In this case 
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For uj m < uj < uj s the relevant expressions are given in eq.( 3.42 ) and ( 3.45 ). 
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For uj < u> m , the relevant expressions are given in eq.(3.60) and (3.69). 
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At a; = u)f-ji we may check, using eq.(l.S) that the spectrum is continuous, up to numerical factors of order 1. 

If lj s > M, then for the range M < uj < uj s the spectrum is approximately the same as the massless spectrum, 
eq.(p~5). 

The determination of the leading contribution is more complicated now, because the ratio of the two terms in each 
range is different and involves also the ratio M-. We postpone such analysis to a more detailed study of specific cases. 

Eqs.(p~15l), ( pU6D may be presented in a more symm etric, perhaps related to some symmetry of the underlying 
physics. We present here only the symmetric form of eq.(4.16). 
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where Mill and Bf 1 ! are numerical coefficients of order one which could be read off eg. (4. 16) 



The spectra of massless and massive particles can be translated into each other by using the relations M/M s = 
(k re /k s ) 2 Zg 2 for the massive case and (k re /k s ) 2 Zg 2 — [ui / ui s ) 2 Zg 2 for the massless case. The substitutions 
(ui/qj s ) 2 z^ 2 <-* M/M s in addition to the substitution accounting for the difference in the energy y/M 2 + up /ui s «-* uj juj s 
translate the results from the massive case into those of the massless case and vice versa. 



V. CONCLUSIONS 



Spectra of produced particles in string cosmology models come in different shapes, as summarized in section 4. 
Even though the background curvature is increasing during the inflationary dilaton-driven phase, spectra are not 
necessarily increasing. This peculiarity comes about because different types of particles couple differently to the 
background curvature and dilaton, and in addition, spectra of massless and massive particles are different. 

During the string phase, the use of the lowest order eguations for the background or for perturbations is guestionable, 
and therefore it might have been possible to doubt the validity and accuracy of the calculations of particle production. 
We have argued that to estimate particle production during the dilaton-driven phase it is not necessary to know the 
details of the evolution during the string phase. If the principle of causality is used, expressed in practical terms by 
the "freezing" of the perturbation amplitude and it's conjugate momentum, it is possible to do reliable calculations. 
Our results provide several explicit checks of this principle, by comparison with other calculations, and therefore give 
further credibility to existing spectral calculations. 

We have not used any explicit symmetry considerations in our calculations. However, it is clear that spectral indices 
and perhaps other properties of the spectra could be analized using an underlying symmetry as in |isf| . 

Among the spectra we found there are some with weak dependence on frequency and some decreasing spectra which 
contain substantial power at large wavelength. Because these spectra depend differently on the basic parameters of the 
models, requiring that they be compatible with astrophysical and cosmological bounds is likely to narrow the allowed 
range of the parameters. Their existence also suggests an obvious source for the observed large scale anisotropy. 

Relic massive particles are an appealing source for cold or hot dark matter. However, to decide on this issue, 
further analysis and input are necessary. To compute relic abundances of particles we need to know more about their 
interactions and also about the late time evolution of the Universe. We are completing such analysis for axions, and 
hope to perform similar analyses for other particles. 

As already mentioned several times, there are obvious improvements and generalizations to our calculations. 
Straightforward improvements, that cannot however be performed in a general way, concern the late time back- 
ground evolution, and the inclusion of effects of particle decay and interactions. Among possible generalizations, we 
view considering other background evolutions as particularly interesting. This was partially done in ]i^ , |l9| and should 
provide information about the robustness of the spectra. 
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